Linearisation instability of gravity waves? 
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Gravity waves in irrotational dust spacetimes are characterised by nonzero magnetic Weyl tensor 
H a i,. In the linearised theory, the divergence of H a i, is set to zero. Recently Lesame et al. [Phys. 
Rev. D 53, 738 (1996)] presented an argument to show that, in the exact nonlinear theory, div H — 
forces H a b = 0, thus implying a linearisation instability for gravity waves interacting with matter. 
However a sign error in the equations invalidates their conclusion. Bianchi type V spacetimes are 
shown to include examples with div H — 7^ H a b- An improved covariant formalism is used to show 
that in a generic irrotational dust spacetime, the covariant constraint equations are preserved under 
evolution. It is shown elsewhere that div H = does not generate further conditions. 

G\ 

C/5 ! I. INTRODUCTION 

O , 

Irrotational dust spacetimes have been widely studied, in particular as models for the late universe, and as arenas 
for the evolution of density perturbations and gravity wave perturbations. In linearised theory, i.e. where the 
irrotational dust spacetime is close to a Friedmann-Robertson- Walker dust spacetime, gravity wave perturbations are 
usually characterised by transverse traceless tensor modes. In terms of the covariant and gauge-invariant perturbation 
formalism initiated by Hawking [Q and developed by Ellis and Bruni these perturbations are described by the 
electric and magnetic Weyl tensors, given respectively by 
0\ 1 , , , , 

O E ab = CacbdU C u a , H ab = ^adetfC^bfU 1 (1) 

\o . 

ON ' where C a bcd is the Weyl tensor, n a bcd is the spacetime permutation tensor, and u a is the dust four-velocity. In the 
so-called 'silent universe' case H a b = 0, no information is exchanged between neighbouring particles, also in the exact 
O , ; nonlinear case. Gravity wave perturbations require nonzero H ab , which is divergence-free in the linearised case ||, 

§, I- 

A crucial question for the analysis of gravity waves interacting with matter is whether the properties of the linearised 
perturbations are in line with those of the exact nonlinear theory. Lesame et al. j| used the covariant formalism and 
then specialised to a shear tetrad, in order to study this question. They concluded that in the nonlinear case, the only 
solutions with div H — are those with H a b = — thus indicating a linearisation instability, with potentially serious 
. £h ! implications for standard analyses of gravity waves, as pointed out in flsl, j ffl. It is shown here that the argument of 
H does not in fact prove that div H — implies H ab = 0. The error in is traced to an incorrect sign in the Weyl 
H [ tensor decomposition (see below) 

The same covariant formalism is used here, but with modifications that lead to simplification and greater clarity. 
This improved covariant formalism renders the equations more transparent, and together with the new identities 
derived via the formalism, it facilitates a fully covariant analysis, not requiring lengthy tetrad calculations such as 
those used in ||. The improved formalism is presented in Section II, and the identities that are crucial for covariant 
analysis are given in the appendix. 

In Section III, a covariant derivation is given to show that in the generic case of irrotational dust spacetimes, the 
constraint equations are preserved under evolution. A by-product of the argument is the identification of the error in 
In a companion paper ||, we use the covariant formalism of Section III to show that when div H = 0, no further 
conditions are generated. In particular, H a b is not forced to vanish, and there is not a linearisation instability. A 
specific example is presented in Section IV, where it is shown that Bianchi type V spacetimes include cases in which 
divH = but H ab ^ 0. 



o 



1 The authors of |3J are in agreement about the error and its implication (private communication). 
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II. THE COVARIANT FORMALISM FOR PROPAGATION AND CONSTRAINT EQUATIONS 

The notation and conventions are based on those of §, §; in particular 8ttG = 1 = c, round brackets enclosing 
indices denote symmetrisation and square brackets denote anti-symmetrisation. Curvature tensor conventions are 
given in the appendix. 

Considerable simplification and streamlining results from the following definitions: the projected permutation tensor 
(compare @, 0), 

Sabc = VabcdU d (2) 

the projected, symmetric and trace-free part of a tensor, 

S<ab> = h a c hb d S( c d) — \S c dh cd h ab (3) 

where h ab = g ab + u a u b is the spatial projector and g a b is the metric, the projected spatial covariant derivative 
(compare @,'|, @), 

B a S c - d e ...f = //.,'//';. • • • h d q h e r ■ ■ ■h f s V b S p -' 1 r ... s (4) 
and the covariant spatial curl of a tensor, 

cur\S ab = e cd{a D c S b) d (5) 

Note that 

S a b = S( ab j curl S ab = curl S <ab> 

since curl ( fh ab ) = for any /. The covariant spatial divergence of S ab is 

(div5) Q = T) b S a b 

The covariant spatial curl of a vector is 

curl5„ =s abc D b S c 

Covariant analysis of propagation and constraint equations involves frequent use of a number of algebraic and 
differential identities governing the above quantities. In particular, one requires commutation rules for spatial and 
time derivatives. The necessary identities are collected for convenience in the appendix, which includes a simplification 
of known results and a number of new results. 

The Einstein, Ricci and Bianchi equations may be covariantly split into propagation and constraint equations ||. 
The propagation equations given in Q for irrotational dust are simplified by the present notation, and become 



p + <dp = (6) 

© + |e 2 = -\p-<J ab o ah (7) 

(Tab + |0O" Qb + (J c<a (T b> C — —E ab (8) 

E a b + QE ab - 3a c<a E b> c = curl H ab - \po ab (9) 

H ab + eH ab - 3a c<a H b> c = -curl^ afc (10) 
while the constraint equations become 

D h cr afc = |D a 6 (11) 

curlCT ah = H ab (12) 

n b E ab = lD a p + e abc a b d H cd (13) 

B b H ab = -e a bc(T b dE cd (14) 



A dot denotes a covariant derivative along u a , p is the dust energy density, & its rate of expansion, and a ab its 
shear. Equations (ph, (|l0|), ( [l3| ) and ( p^ ) display the analogy with Maxwell's theory. The FRW case is covariantly 
characterised by 

T) a p = = D a 9 , Oab = Eab = H ab = 



2 



and in the linearised case of an almost FRW spacetime, these gradients and tensors are first order of smallness. 

The dynamical fields in these equations are the scalars p and O, and the tensors o~ a b, E a b and H a b, which all satisfy 
S a b = S <a b>- The metric h a b of the spatial surfaces orthogonal to u a is implicitly also involved in the equations as a 
dynamical field. Its propagation equation is simply the identity h a b = 0, and its constraint equation is the identity 
D a /if, c = - see ( A5). The Gauss-Codacci equations for the Ricci curvature of the spatial surfaces M 



\R*h a b 
R* 



-|e 2 



a ab <J ab + 2p 



have not been included, since the curvature is algebraically determined by the other fields, as follows from (f 



E„ 



0~ ca o~b 



(15) 



(16) 



The contracted Bianchi identities for the 3-surfaces 



D b R* ab = lD a JT 



reduce to the Bianchi constraint ( J13| ) on using (jil|), (12) and the identity (A15) in ( |l5| ) and (|16|). Consequently, these 
identities do not impose any new constraints. 

By the constraint (112 V one can in principle eliminate H a b. However, this leads to second-order derivatives in the 
propagation equations (S) and (|Io|). It seems preferable to maintain H a b as a basic field. 

One interesting use of (n~2) is in decoupling the shear from the Weyl tensor. Taking the time derivative of the 
s hear propagation equation^) , using the propagation equation (0) and the constraint (|l2) , together with the identity 
( Aisp , one gets 



-D 2 a a6 



0~ab 



TQa a b — hQ(T ab 



46cr c<a (Tb> c + Q(T cd a c<a (Jb>d - 2& ae <7deh c <a&b> c + 4cr c<a (Tb> 



|D <a D c crb >c 

de. 



(17) 



where D 2 = D a D a is the covariant Laplacian. This is the exact nonlinear generalisation of the linearised wave equation 
for shear perturbations derived in In the linearised case, the right hand side of ([l7j ) vanishes, leading to a wave 
equation governing the propagation of shear perturbations in an almost FRW dust spacetime: 



-D 2 a ab 



0~ab 



D <a DV 6>c w 



As suggested by comparison of ( J12| ) and (J14|), and confirmed by the identity (A17), div curl is not zero, unlike its 
Euclidean vector counterpart. Indeed, the divergence of ( |l2"| ) reproduces (fu]), on using the (vector) curl of ( pi] ) and 
the identities (DO), dAS) and (|A17|): 



div (12) and curl (|ll|) -> (|lj) 

Further differential relations amongst the propagation and constraint equations are 

curl (8) and ([□]) and (|l|) and (||)' -> © 
grad (|) and div (|) and (|ll|) and (|ll|)' and (|l|) -» @ 



(18) 

(19) 
(20) 



where the identities (fV|), ( |A13[ ), ( |A15D , ( |A16| ) and ( |A18| ) have been used. 

Consistency conditions may arise to preserve the constraint equations under propagation along u a ||, jij. In the 
general case, i.e. without imposing any assumptions about H a b or other quantities, the constraints are preserved 
under evolution. This is shown in the next section, and forms the basis for analysing special cases, such as diviJ = 0. 



III. EVOLVING THE CONSTRAINTS: GENERAL CASE 

Denote the constraint equations ([ll]) — ( |l4| ) by C A = 0, where 

C A = (DV Qb - |D a 8 , curiae - H ab , ■ ■ •) 

and A — 1, • ■ • , 4. The evolution of C A along u a leads to a system of equations C A = J rA (C B ), where T A do not 
contain time derivatives, since these are eliminated via the propagation equations and suitable identities. Explicitly, 
one obtains after lengthy calculations the following: 
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C 1 a — — 6C 1 a + 2e a bc ab d C 2 c d — C 3 a 

C 2 ab = — @C 2 a f, — £ C!i ( a Crh) c C 1 d 

r 3 — —iPir 3 4- x n h r 3 — Lnr 1 

^ a 3 yJL' a i 2 a ^ 2 ^ a 



c 4 



I- ^EaC^b — e a bc Et, d C 2 c d + ^curlC 4 a 
l~ f H b C X b — £a bc Hb d C 2 c d — 5CurlC 3 a 



(21) 
(22) 

(23) 

(24) 



For completeness, the following list of equations used in the derivation is given: 

Equation requires @, (|A^)jf7|), (§), @, @, (|l|), ( |A15| ) - where ( |A1~5| ) is needed to eliminate the following 
term from the right hand side of (pi|): 

s abc (7 b dCurla cd - a bc T> a a b c 
+ a bc B c a ab + \oacV h o bc = 



Equation (^) requires ( A18 ), (||), (|ic|), (|Tl|), (|l2|), (A4) - where (A4) is needed to eliminate the following term from 
the right hand side of (|22j): 

s cd (a {D c [a b fa d e ] + D e [a b) d a c e ] } = 



Equation @ requires (^T|), (§, (|), ©, flATp) , Q, ([Tl]), (|l3|), Q, ( |AT5| ) - where (^T|) is needed to eliminate 
the following term from the right hand side of (|23|): 

\a ab D c E bc + e abc E b d curl a cd 

+ e abc a b d curl E cd + \E ab D c a bc + E bc D b a ac 

+ a bc D b E ac - D a (a bc E bc ) = 



Equation @ requires (§, (|), ©, @), gl|L_P), @, @, ©■ 

In ||] , a sign error in the Weyl tensor decomposition ( |A€ ) led to spurious consistency conditions arising from the 
evolution of ([H]), fll^)- The evolution of the Bianchi constraints (|l3|), ( p"4| ) was not considered in ||. 

Now suppose that the constraints are satisfied on an initial spatial surface {t — to}, i.e. 



= 



(25) 



where t is proper time along the dust worldlines. Then by ( |21| ) - (24), it follows that the constraints are satisfied for 
all time, since C A = is a solution for the given initial data. Since the system is linear, this solution is unique. 

This establishes that the constraint equations are preserved under evolution. However, it does not prove existence 
of solutions to the constraints in the generic case — only that if solutions exist, then they evolve consistently. The 
question of existence is currently under investigation. One would like to show explicitly how a metric is constructed 
from given initial data in the covariant formalism. This involves in particular considering whether the constraints 
generate new constraints, i.e. whether they are integrable as they stand, or whether there are implicit integrability 
conditions. The relation ( |l8| ) is part of the answer to this question, in that it shows how, within any {t = const} 
surface, the constraint C 4 is satisfied if C 1 and C 2 are satisfied. Specifically, (Uq) shows that 



C 4 a = icurlC^ - D b C 2 Q0 



(26) 



Hence, if one takes C 1 as determining gradG, C 2 as defining H and C 3 as determining grad/?, the constraint equations 
are consistent with each other because C 4 then follows. Thus if there exists a solution to the constraints on {t = to}, 
then it is consistent and it evolves consistently. 

In the next section, Bianchi type V spacetimes are shown to provide a concrete example of existence and consistency 
in the case 



divE / ^ cur\E, divH = ^ curliJ , grad p = = grad 9 
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IV. SPACETIMES WITH DIVtf = / H 



Suppose now that the magnetic Weyl tensor is divergence-free, a necessary condition for gravity waves: 

divtf = [a,E}=0 (27) 

where [S, V] is the index-free notation for the covariant commutator of tensors [see (|A2|)], and the equivalence follows 
from the constraint (|l4|). Using the covariant formalism of Section III, it can be shown Q that ( |27j ) is preserved 
under evolution without generating further conditions. In particular, ( p7j ) does not force H ab = - as shown by the 
following explicit example. 

First note that by and (||): 

R* ab = \R*h ab [a,R*}=0 div H = 

i.e., irrotational dust spacetimes have divH = if R* b is isotropic. 

Now the example arises from the class of irrotational spatially homogeneous spacetimes, comprehensively analysed 
and classified by Ellis and MacCallum According to Theorem 7.1 of |L3]], the only non-FRW spatially homo- 

geneous spacetimes with R* b isotropic arc Bianchi type I and (non-axisymmctric) Bianchi type V. The former have 
H ab = 0. For the latter, using the shear eigenframe {e a } of |]l3| 

a a& = (7 22 diag(0, 0,1,-1) (28) 

Using ^) and @ with @, one obtains 

E a b = j9^tt6 — Oc<aO~b> C 

= \(J 2 1 diag (0, 2(722 , 6 - (722, -6 - (722) (29) 

in agreement with p3|.p| 

The tetrad forms of div and curl for type V are (compare |M| ): 

D b 5 afc = a 6 5 Q b - 3a b 5 a6 (30) 
curlS'afc = e cd ( a d c S b - } d + e cd ( a S b ) c a d (31) 

where S^f, = S <ab> , a b = aSi, 1 (a is the type V Lie algebra parameter) and d a f is the directional derivative of / along 
e a . Using p|) and A 



H a b = curl a a 

7220( a O b ) 



2aa 22 6 (a 2 6 b) 3 (32) 



Hence: 

Irrotational Bianchi V dust spacetimes in general satisfy div H — ^ H. 



Using (|28|) — d32|), one obtains 

B b H ab = (33) 

curliJ afc = -a 2 cr ab (34) 

curl curl H ab = —a 2 H ab (35) 

D b ^ afc = -a bc a bc a a (36) 

curl£ a6 = i9i/ ab (37) 



Although ( |33| ) is a necessary condition for gravity waves, it is not sufficient, and (g9|) and (|32| ) show that E ab and 
H ab decay with the shear, so that the type V solutions cannot be interpreted as gravity waves. Nevertheless, these 
solutions do establish the existence of spacetimes with diviJ = ^ H. 

This supplements the known result that the only spatially homogeneous irrotational dust spacetimes with H a b = 



are FRW, Bianchi types I and VI_i (n a a = 0), and Kantowski-Sachs p5|. When H ab = 0, (E8h and (32) show that 



2 Note that E a b in 113] is the negative of E a b defined in (hi). 
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<Jab = 0, in which case the type V solution reduces to FRW. 

A final remark concerns the special case H ao = 0, i.e. the silent universes. The considerations of this paper show 
that the consistency analysis of silent universes undertaken in fl(| needs to be re-examined. This is a further topic 
currently under investigation. It seems likely that the silent universes, in the full nonlinear theory, arc not in general 
consistent. 
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APPENDIX: COVARIANT IDENTITIES 



In this appendix 

S a b = S <a b> and V a b — V <a b> 
Using the properties of r] a bcd (see ||) and h a b, one can derive 

e abc e def = 3\h^ d h b e h c ^ (Al) 

SabcS b d V cd = \Eabc [S, Vf C (A2) 

where [S, V] is the covariant commutator of rank- 2 tensors. Using a tetrad that diagonalises S ao , one may prove the 
further covariant identities: 

e a bcS h P S%V c « = -Sa b e bcd S c p V dp (A3) 
curl (S 2 ) ab = e cd(a B e {S b) c S d e } (A4) 

where (S ) a b = S a c S c b- Furthermore 

V a h bc = = B a e bcd , h ab = Q = e a bc (A5) 

where the time derivative identities depend on u a — 0, while the spatial derivative identities do not. 
The curvature tensor is 

Rabcd = Cabcd + g a [cRd]b — 9b[ C Rd]a ~ \R9abcd 

where g a bcd = 9ac9bd - 9ad9bc, and 

Cabcd — {.9abpq9cdrs ^JabpqVcdrs)'^'^' 

~ {riabpqgcdrs + 9abp q r]cdrs)u P U r H QS (A6) 

(see 0, corrected here and in agreement with &). By the field equations for dust, the Ricci tensor is 

Rab = R C acb = \p{u a Ub + h ab ) 



and, together with (A6), this gives 

Rabcd — 2 (h a [ c + u a U[ c ) E d ] b + 2E a [ c (h d ^ b + u d ]U b ) 

+ 2£abeU[ c H d f + 2e cde U[ a H b f 

+ \p (h a[c u d] u a - h b[c u d] u a + 2h a[c h d]b ) (A7) 

The Ricci identities f- : [ a b] — and 

2^0; [be] = R d abcY d 



G 



2W, 



ab\ [cd] 



R e acdW t 



R e bcdW a 



together with (A7) and 



u a -b 



\Qh ab + <J ab , U a = = UJab 



lead to the following crucial identities: 



(Da/)' 


= D a /- 


D [a D b] / 


= 


(D a s b y 


= B a S b - 




= (i° 2 - 




+ s [b \ 


(VaSbcY 


= D a 5b c 


(v b s ab )' 


= D b S ab 


D [a D 6] S cd 


= 2 (ie 2 



i9D a / - a a b l 



^•au a Dbc — CT a UdDbc -t 6n a bde{b^c) 

D b 5 , ab — \&D b S ab — o~ bc D c S a b + e a b c H b dS cd 



2S { J C [E b] d) - |6<7 6] d )} 



2^a (C {^] e -| Q(7 %}^ )e 



E cd ) 



Then ( Ra )— ( A14) imply the further important identities 

e abc S b p cm\V cp = 2S bc D [a V b]c - \S ab D c V bc 
curl (fS ab ) = f curl S ab + e cd(a S b) d V c f 
D b curlS Qb = ±e Qbc D b (V d S cd ) + e abc S b d (±6a cd 

— <T ab e (T ce b d 
(curlS'ah) = cm\S a b - ■|9curlS' ob 

— ( J e C £cd(a^ >e Sb) d + 3H c<a S b > C 

curl curl S ab = -B 2 S ab + §D <Q D c S b>c + (p- ±6 2 ) S at 
+ 3S c<a {E b> c — |Ocr b> c } + a cd S c <Jab 

— S cd a ca (Jbd + cr cd (J c ( a Sb)d 



D 2 (Da/) =D a (D 2 /)+i?:"D b / 



(A8) 
(A9) 
(A10) 

(All) 
(A12) 

(A13) 



(A14) 

(A15) 
(A16) 

(A17) 

(A18) 



(A19) 
(A20) 



where D 2 = D a D a is the covariant Laplacian. 
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